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Recent work on QCD-like theories has shown that the addition of adjoint fermions obeying periodic 
boundary conditions to gauge theories on R 3 x S 1 can lead to a restoration of center symmetry and 
confinement for sufficiently small circumference L of S 1 . At small L, perturbation theory may be 
used reliably to compute the effective potential for the Polyakov loop P in the compact direction. 
Periodic adjoint fermions act in opposition to the gauge fields, which by themselves would lead to a 
deconfined phase at small L. In order for the fermionic effects to dominate gauge field effects in the 
effective potential, the fermion mass must be sufficiently small. This indicates that chiral symmetry 
breaking effects are potentially important. We develop a Polyakov- Nambu-Jona Lasinio (PNJL) 
model which combines the known perturbative behavior of adjoint QCD models at small L with 
chiral symmetry breaking effects to produce an effective potential for the Polyakov loop P and the 
chiral order parameter •tptp. A rich phase structure emerges from the effective potential. Our results 
are consistent with the recent lattice simulations of Cossu and D'Elia, which found no evidence for a 
direct connection between the small-L and large-L confining regions. Nevertheless, the two confined 
regions are connected indirectly if an extended field theory model with an irrelevant four-fermion 
interaction is considered. Thus the small-L and large-L regions are part of a single confined phase. 



I. INTRODUCTION 



Recent progress in the study of QCD-like gauge theories has revealed that a confined phase can exist under certain 
conditions when one or more spatial directions are compactified and small [T] [2]. This is surprising, because a 
small compact direction in Euclidean time gives rise to a deconfined phase for SU(N) gauge theories. It is also 
intriguing, because one or more small compact directions give rise to a small effective coupling constant if the theory 
is asymptotically free. Thus we now have four-dimensional field theories in which confinement holds, and holds under 
circumstances where semiclassical methods may be reliably applied. In addition, the existence of confinement with one 
or more compact directions suggests that these new theories may make possible the construction of large- N models 
in small space-time models, finally realizing the potential of the Eguchi-Kawai large- TV reduction |3] SJ [5] . 

When one or more directions are compact, the perturbative contribution of the gauge fields to the effective potential 
favors the deconfined phase. This leads to the well-known deconfinement transition in finite-temperature gauge 
theories. Thus it is necessary to modify the gauge theory in some way to obtain confinement with small compact 
directions. At present, there are two methods known for achieving this. The first method directly modifies the gauge 
action with terms non-local in the compact direction(s) [2], while the second adds adjoint fermions with periodic 
boundary conditions in the compact direction(s) [T], which is our subject here. 

The term adjoint QCD is often used to refer to vector gauge theories based on the gauge group SU(N) where, in 
addition to the gauge fields, fermions in the adjoint representation are fundamental fields. We will use Nf to denote 
the number of flavors of adjoint Dirac fermions. With this notation Nf — 1/2 corresponds to TV = 1 supersymmetry 
when the mass of the fermion is taken to zero. If Nf < 2, the theory will be asymptotically free. Here we will consider 
models on i? 3 x S 1 that give rise to a confining phase if the circumference L of S 1 is sufficiently small and the mass 
m of the adjoint fermions is sufficiently light. 

Confinement in SU(N) gauge theories is associated with an unbroken global center symmetry, which is Z(N) for 
SU(N). The deconfined phase is associated with the breaking of this Z(N) symmetry, which occurs spontaneously 
at high temperatures [6, [7]. The order parameter for Z(N) breaking in the compact direction is the Polyakov loop 
P, which is the path-ordered exponential of the gauge field in the compact direction 



P (x) = V exp 



i / dx^A^ (x) 
Jo 



(1) 
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The trace of P in a representation R represents the insertion of a heavy fermion in that representation into the 
system. Unbroken Z(N) symmetry implies (TrpP) = in the confined phase, and correspondingly (TrpP) ^ holds 
in the deconfined phase where Z(N) symmetry is broken. The character Tr^P of each irreducible representation 
of SU(N) transforms as TtrP — > z fe Tr^P under P — > zP for some k <E {0, ..,N — 1}. In the confined phase, the 
expectation value (Tr^jP) is for all representations that transform non-trivially under Z(N), i.e., have k ^ 0. In 
the case where bosons have periodic boundary conditions and fermions have antiperiodic boundary conditions in the 
compact direction, the transfer matrix in that direction is positive-definite and the system can be interpreted as being 
at temperature T = = L~ l . In gauge models where all physical fields have zero 7V-ality, the free energy F for a 
single test fermion in the fundamental representation is related to the expectation value of Polyakov loop 

e -PF =< TrpP > _ ( 2 ) 

The restoration of Z(N) symmetry for small L and small m when the adjoint fermions have periodic boundary 
conditions is seen from the behavior of the one-loop effective potential. This calculation is a variant of the calculation 
of the effective potential at finite temperature j6] |7] |8]. The sum of the effective potential for the fermions plus that 
of the gauge bosons gives 



Vi-ioop (P, L, m, N f ) = -2=4 



Tr A P n 



L 4 

71 = 1 



2N f L 2 m 2 K 2 (nLm) - \ 



(3) 



Note that the first term in brackets, due to the fermions, is positive for every value of n, while the second term, due 
to the gauge bosons, is negative. In the limit m — ► 0, this becomes 

v^ loop * Ef,^-^) 2 -^^ 1 [4>i - M 2 (4>i - <fe - 2tt) 2 ] (4) 

where we have written P as Pjk = 5jke l ^ j in a gauge where the background field A4 is diagonal and independent of 
X4. If mL is sufficiently small, this effective potential has a global minimum when the Polyakov loop eigenvalues are 
uniformly spaced around the unit circle. This is the unique Z(iV)-symmetric solution for P. On the other hand, if 
mL is sufficiently large, the fermion contribution to Vi-i 00 p ls negligible compared to the gauge boson contribution, 
and the system will be in the deconfined phase. Interestingly, there can be intermediate phases between the confined 
and deconfined phases as m is varied for N > 3 019]. It is also interesting to note that the one-loop effective potential 
identically vanishes when Nf = 1/2 in the m = limit, consistent with the supersymmetry of the model |10l [TT| . 

It is clearly crucial that the fermion mass m be sufficiently small in order for confinement to be restored at small L. 
Experience with phenomenological models |12[ fT3] suggests that in fact it is the constituent mass which is relevant in 
determining the size of the fermionic contribution to the effective potential for P. The potential importance of chiral 
symmetry breaking and restoration is underscored by lattice results for adjoint SU(3) fermions at finite temperature, 
where the fermions are antiperiodic in the Euclidean timelike direction. In this case, there is a large separation between 
the deconfinement temperature and the chiral symmetry restoration temperature, with T c /Td ~ 7.8(2) |14l Il5] . 

In order to explore the interrelationship of confinement and chiral symmetry breaking, we use a generalization 
of Nambu-Jona Lasinio models known as Polyakov-Nambu-Jona Lasinio (PNJL) models [13]. In NJL models, a 
four-fermion interaction induces chiral symmetry breaking. There has been a great deal of work on NJL models, 
both as phenomenological models for hadrons and as effective theories of QCD |16| IT7] . NJL models have been 
used to study hadronic physics at finite temperature, but they include only chiral symmetry restoration, and do 
not model deconfinement. This omission is rectified by the PNJL models, which include both chiral restoration and 
deconfinement. The earliest model of this type was derived from strong-coupling lattice gauge theory [12j, but later 
work on continuum models have proven to be extremely powerful in describing the finite-temperature QCD phase 
transition [13]. In PNJL models, fermions with NJL couplings move in a non-trivial Polyakov loop background, and 
the effects of gluons at finite temperature is modeled in a semi-phenomenological way. We will develop a model of 
this type for both fundamental and adjoint fermions below. 

These models will show why mL, with m interpreted as a constituent mass, is the key parameter. Because we 
expect m at worst to stay constant as L is decreased, the condition for confinement that mL is small can always be 
met at sufficiently small L. It is not a priori obvious what sufficiently small means, given the persistence of chiral 
symmetry breaking in the deconfined phase. Recent lattice simulations by Cossu and D'Elia [183 nave confirmed the 
existence of the small-L confined region in SU(3) lattice gauge theory with two flavors of adjoint fermions, and we will 
focus on this case in our analysis. Even if the small- L confined region exists and is accessible in lattice simulations, it 
is not necessarily the same phase as found for large L. Put slightly differently, we would like to know if the small- L 
and large-L confined regions are smoothly connected, and thus represent the same phase. Our main result will be 
a phase diagram for adjoint periodic QCD for all values of L, obtained using a PNJL model. On the way to this 
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goal, we will use as tests of our model both standard QCD with fundamental fermions and adjoint QCD with the 
usual antiperiodic boundary conditions for fermions. Our principal tool will be the effective potential for the chiral 
symmetry order parameter ipip and the deconfinement order parameter P. 

In section II the fermionic contribution to the effective potential is derived, and section III discusses the gluonic 
contribution. The complete effective potential is tested with fundamental fermions in section IV. Sections V and VI 
discuss adjoint fermions with antiperiodic and periodic boundary conditions, respectively. A final section summarizes 
our conclusions. 



II. FERMIONIC CONTRIBUTION TO EFFECTIVE POTENTIAL 



We start our calculations by constructing the contributions to the effective potential from both fundamental and 
adjoint fermions. NJL models use purely fermionic interactions as a proxy for the gauge theory interactions that give 
rise to chiral symmetry breaking. Typically, the relation between the gauge theory and an associated NJL model 
is fixed by matching important hadronic parameters such as f n . In the case of NJL models at finite temperature, 
and particularly PNJL models, it is common to assume that the NJL model parameters are fixed by T = hadronic 
parameters and remain constant as T is increased, at least up to the deconfinement temperature. At substantially 
higher temperatures, it is possible that the parameters of the NJL model might have to be adjusted to reproduce 
properties of the gauge theory at those higher temperatures. This possibility is particularly acute in adjoint QCD, 
where the scale of chiral symmetry restoration is almost a factor of ten higher than the deconfinement scale |14l [T5] . 
Therefore, it is desirable to consider the phase diagram of the appropriate PNJL model over a range of parameters, 
holding open the possibility that the PNJL parameters are varied with (3 or L. 

We take the fermionic part of the Lagrangian of our PNJL model to be |13l [TBI IT?] 

L F = ^(< 7 . D - to ) V + y [(?M» 2 + (^7sA» 2 ] + g D [det V; (1 - 7s) V> + h.c] (5) 

where ip is associated with Nt flavors of Dirac fermions in the fundamental or adjoint representation of the gauge 
group SU(N). The A a 's are the generators of the flavor symmetry group U(Nf); g$ represents the strength of the 
four-fermion scalar-pseudoscalar coupling and gp> fixes the strength of an anomaly-induced term. For simplicity, we 
take the Lagrangian mass matrix too to be diagonal: (mo)j k — mojSjk- If too and gp, are taken to be zero, Lp is 
invariant under the global symmetry U(Nj)p x U(Nf)n. With a non-zero, flavor-independent mass to and gp ^ 0, 
the symmetry is reduced to SU(Nf)v x E7(l)v- The covariant derivative _D M couples the fermions to a background 
Polyakov loop via the component of the gauge field in the compact direction. 
We will use the identity 

AgfAg, = 25 a 5 jh (6) 
for the generators of U(Nf) to rewrite the scalar-scalar four-fermion coupling as 

55 (^V-i)] • (7) 

The terms with i — j 



as 



E(^) 2 (8) 



give rise to an interaction that does not mix flavors. These terms are of order jV 2 , as opposed to the terms with 
i j, which are only of order N. Keeping only the terms with i = j is the Hartree approximation, but the i ^ j 
contribution to the one-loop vacuum energy is canceled by the pseudoscalar interaction in any case. We will make a 
similar approximation for the gp, interaction. 
The partition function associated with Lp is 

Z F = j [#] [<ty]e* (9) 

and depends implicitly on the background gauge field via the covariant derivative. We introduce a set of auxiliary 
scalar fields, <jj and a„ for each flavor j, and write the partition function as 



Z = J [dif] [dip] [da] [da] exp i J d^x I Lp + otj (aj — ipjipj 



(10) 
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where the <x,'s provide functional 5 functions that set aj equal to ipjipj. This allows us to rewrite Z in the form 



Z = j [dtp] [dtp] [da] [da] exp J i J d 4 



i>i (h-D- m 0j - aj ) + 9stf + 2g D JJ trj + a j°j 



We proceed by integrating over the fermion fields to obtain 



j [da] [da] 



Z= I [da] [da] e lS ^f 



where the effective action S e ff is given by 



S e ff = —i In [det (ry • D — m j — aj)] + / d 4 x 



+ 2 9d Yl a i ■ + /Z a i a i 



■ (11) 



(12) 



(13) 



We look for the stationary saddle points of the effective action regarded as a function of the atj's. The expression to 
be minimized has the form 



zZ F ^ + zZ a i a o ~V(a) 



(14) 



whose saddle point solutions satisfy 



9Fi 



(15) 



dV n 

a >-do-r {) 



(16) 



which reduces to 



This in turn is equivalent to extremizing 



of) 



*j = 0. 



3 J 



(17) 



(18) 



In this particular case, S e ff is reduced to 



~>eff 



+ / d 4 x 



dot | 17 • D - m 0j + 2g s uj + 2g D ]J a k 



We see immediately that each flavor behaves as if it has a mass given by rrij — m j — 2g$aj — 2gDTik^j a k- The 
boundary conditions in the compact direction enter the effective potential through the fermion determinant, so it is 
through the mass mj, the constituent mass, that the Polyakov loop and the chiral order parameters aj are coupled. 
It is generally convenient to use the language of finite temperature to describe both the case of finite temperature, 
= T > 0, with antiperiodic boundary conditions, and the case of a periodic spatial direction, L < oo. The 
fermionic contribution to the effective potential has both a zero-temperature and a finite-temperature contribution. 
The zero-temperature part consists of a potential term, given by 



J2gs<?- + 2g D (N f -l)l[a, 



(19) 
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as well as a contribution from the fermion functional determinant. This is formally given by 



where oJf! — i /fc 2 + to 2 , and the constant du is the dimensionality of the color representation, N for the fundamental 



and N 2 — 1 for the adjoint. This contribution, representing a sum of one-loop diagrams, is divergent, and requires 
regularization. Although many schemes have been used [16], a non-covariant three-dimensional cutoff has been used 
most often. With this regularization, the zero-temperature part of the fermion determinant is given by 

N f 2A y^A 2 + to 2 (2A 2 + to 2 ) + m) [log (to 2 ) - 2 log (a + ^/A 2 + to 2 )] 
-2^E 3^2 (21) 

where A is the momentum space cutoff. Thus the total T — fermionic contribution to the effective potential is 
V P o(m,m ,) = E ffsfT i + 2g D (N f - 1) JJctj 
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w / 2A, /A 2 + to 2 (2 A 2 + m 2 ) + m 4 [log (m 2 ) - 2 log (a 



32tt 2 



(22) 



In PNJL models, the finite-temperature contribution from the fermion determinant depends on the background 
Polyakov loop. It is convenient to work in a gauge where the temporal component of the background gauge field, 
A±(£, t), is constant and diagonal. The covariant derivative then becomes 7 • D = 7 • d — i^A^. The one-loop free 
energy of fermions in a representation R of SU(N) gauge theory with zero chemical potential can be written as 

V FL (P,m) = -2j2Tr R [^ J ^l n (l T Pe~ L ^ ) + h.c] (23) 

where the minus sign is used for periodic boundary conditions and plus for antiperiodic. In the latter case, L is 
replaced by (3. The potential term Vft has a series expansion in terms of modified Bessel functions 

^ ^ "0 = £ Sir E {±irT J RPn K 2 (nL/TTij ) (24) 

j n=l 

which is rapidly convergent for all values of the mass j8]. 

Although NJL models display chiral symmetry breaking and are plausibly related to QCD-like theories, they do not 
capture exactly all the features of the gauge theory to which they are related. There are many ways of associating a 
given NJL model with a gauge theory, and there is freedom in choosing both regularization scheme and parameters of 
the two theories to match. Typically, the hadronic parameters used represent hadronic physics at zero temperature, 
and those parameters are held fixed up to the chiral transition and perhaps beyond. For the case of N = 3 models 
with two and three flavors in the fundamental representation, appropriate choices of Too, gs and A can capably model 
many features of hadron physics. In what follows, we will take Nf = 2, and take the masses moj to be equal to a 
common mass which we also write as mo- In this case, the contribution to S e ff from gs and go has the same form. 
It is convenient to take go = 0, and also to write the common constituent mass as to — m Q — 2gg(T [16, T7j- In PNJL 
models, additional information about confinement at non-zero temperature must be added. As will be discussed 
below, we will use as input the deconfinement temperature Td for the pure gauge theory. 

At non-zero temperature, the observables of a given gauge theory will evolve with T according to the finite- 
temperature renormalization group, but the parameters of an associated NJL or PNJL model do not naturally evolve 
in a related way. If we follow the behavior of both the gauge theory and an associated PNJL model over a large 
scale of temperatures, it may be necessary or desirable to consider the parameters of the PNJL model as varying with 
T. For example, Unsal has developed a comprehensive scenario for periodic adjoint fermions with too = in which 
chiral symmetry is spontaneously broken by monopole-induced anomaly terms [TJ[T9]. If this interaction is sufficiently 
weak, it may not induce chiral symmetry breaking, raising the possibility of confinement without chiral symmetry 
breaking. There is also the possibility of modifying the strength of chiral symmetry breaking by adding additional 
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Figure 1: The normalized order parameter {TrpP)/3 for pure SU(3) gauge theory as a function of temperature. The decon- 
fmement transition is first-order, and occurs at Td = 270 MeV. 



couplings compatible with all symmetries have been added. The easiest construction of such an extended model might 
be obtained by adding to a gauge theory a non-renormalizable four-fermion coupling of exactly the same form as in 
an NJL model, but with opposite sign for the coupling: 



5gs_ 
2 



(25) 



If we imagine that the gauge theory is associated with an NJL model Lnjl with a certain value of g$, then the 
additional term shifts the NJL Lagrangian as 



L 



NJL 



^NJL 



(^A Q V) 2 +(^7 5 A a V) 2 



(26) 



where one expects 8g' s = 5g$ for a sufficiently small perturbation. In principle, a lattice version of the extended gauge 
theory could be simulated, and the results combined with analytical results for the NJL model to determine 5g' s as 
a function of Sgs- In fact, lattice gauge theories with additional non-renormalizable four-fermion terms added have 
already been used in the study of finite temperature gauge theories |20| [21] [22]. For the two-flavor case we discuss, it 
is useful to treat all parameters of the model as potentially varying. 



III. GLUONIC CONTRIBUTION TO EFFECTIVE POTENTIAL 



For gauge bosons, the one-loop finite temperature free energy in a background Polyakov loop is given by an 
expression similar to the one for fermions. The boundary conditions for the gauge bosons are periodic in all cases 
considered here, so L and f3 may be used equivalently in the gluonic sector. We have 



Vg— 1 loop 



1 f r/ 3 k 

(P) = 2Tr A [ I J —ln(l-Pe- L ^)} (27) 



where we have inserted a mass parameter in Clf. — \Jk 2 + M 2 for purely phenomenological reasons explained below. 
If we take the zero mass limit of the series expansion, we get the standard expression for the one-loop gauge boson 
free energy j6] [7] 

M^o n 2 L 2 ^— ' n 2 7r 2 L 4 ^ n 4 
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This infinite series may be summed exactly, giving an expression for V g -ii oop proportional to L~ 4 involving the fourth 
Bernoulli polynomial. If we also retain the next-order term in a high-temperature expansion, proportional to M 2 /L 2 , 
we obtain a useful phenomenological model for the deconfinement transition in pure gauge theories [23j. The potential 
takes the form 

-k a L^ ' n 4 atx L '-^ n z 

n—l n— 1 

We stress that the mass parameter M should not be interpreted as a gauge boson mass, nor do we limit ourselves to 
ML <C 1- The crucial feature of this potential is that for sufficiently large values of the dimensionless parameter ML, 
the potential leads to a Z(iV)-symmetric, confining minimum for P [9] [23]. On the other hand, for small values of 
ML, the pure gauge theory will be in the deconfined phase. It will be important later that V g is a good representation 
of the gauge boson contribution for high temperatures; in other PNJL models, the gauge boson contribution has 
sometimes been chosen so as to be valid over a more narrow range of temperatures. Both of the infinite sums can 
be carried out exactly, giving a closed form for V g as a function of the angles <f>j in terms of the fourth and second 
Bernoulli polynomials |23] . 

In the Polyakov gauge, the Polyakov loop in the fundamental representation of SU(3) can be written as Pjk — 
exp(i(f)j)5jk with two independent angles. With the use of Z(3) symmetry, it is sufficient to consider the case where 
(TtfP) is real. Thus we consider only diagonal, special-unitary matrices with real trace, which may be parametrized 
by taking 4>i = <f>> 4>2 = — 4>i an d 4>3 = 0, or P = diag [e 4( ^, e~ i( ^, l] with < <f> < it. The unique set of Z(3)-invariant 
eigenvalues are obtained for <f> = 2ir/3. For SU(3), V g takes the form: 

v <p\ - f' 3 ^ M ' 2 j_ 1 ^ 1350 4 - 3OO7T0 3 + 180^ 2 - 16^ 

Vg Kn ~ \2n 2 7r 3/ L 2 L 4 \ 90^ 2 ) (3Uj 

We will set the mass scale M by requiring that V g yields the correct deconfinement temperature for the pure gauge 
theory, with a value of T<j w 270 MeV. This gives M = 596 MeV [23 . The pressure p is given by the value of — V g 
at the minimum of the potential. The behavior obtained for the pressure, the energy density e, and the so-called 
interaction measure A = (e — 3p) /T 4 are all consistent with lattice simulations for T > T^. As shown in figure [l] the 
order parameter TrP for the deconfinement transition jumps at Td, indicating a first-order deconfinement transition 
for SU(3). 



IV. FUNDAMENTAL FERMIONS 



As a test of all the components of the effective potential we have assembled, we consider the case of two flavors 
of fundamental fermions at finite temperature. The fermions obey antiperiodic boundary conditions, so we identify 
the compact circumference L with the inverse temperature /3 = 1/T. A very common choice of zero-temperature 
parameters for two degenerate light flavors is m = 5.5MeV, A = 631.4 MeV, and g$ — 2 x 5.496GeV~ 2 [13j [T7] . 
In figure [2] we show the expectation value of the Polyakov loop in the fundamental representation TrpP and the 
constituent mass m as functions of the temperature. The behavior in the crossover region is very similar to the 
results of Fukushima [13], and shows the explanatory power of PNJL models. The constituent mass m is heavy at low 
temperatures, due to chiral symmetry breaking. The larger the constituent mass, the smaller the Z(3) breaking effect 
of the fermions. On the other hand, a small value for (TrpP) reduces the effectiveness of finite-temperature effects 
in restoring chiral symmetry. These synergistic effects combine in the case of fundamental representation fermions to 
give a single crossover temperature at which both order parameters are changing rapidly, in agreement with lattice 
simulations. 



V. ADJOINT FERMIONS WITH ANTIPERIODIC BOUNDARY CONDITIONS 

Adjoint SU(3) fermions at finite temperature show a completely different behavior in lattice simulations from 
fundamental fermions. Because the adjoint fermions respect the Z(3) center symmetry, there is a true deconfinement 
transition where Z(3) spontaneously breaks. Lattice simulations have shown that chiral symmetry is restored at 
a substantially higher temperature than the deconfinement temperature |14l [15]. Unlike the case of fundamental 
fermions, there are no comprehensive T — lattice results for adjoint SU(3) fermions, no established parameter sets 
for NJL models and, of course, no experimental results to provide guidance. We will again consider the case Nt = 2 
with degenerate fermion masses. The T = parameters needed are gs and A. The mass parameter M is again 
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Figure 2: The constituent mass m and (TtfP) for two-flavor QCD with fundamental representation fermions with antiperiodic 
boundary conditions as a function of temperature. The order parameters are normalized by dividing by their values at T — 
and T = oo, respectively. 

fixed by the deconfinement transition in the pure gauge theory, and gp is again absorbed into g$- Rather than work 
directly with g$, we will consider the dimensionless coupling k — gsA 2 . A given ratio of m(T = 0)/A determines the 
value of k, and vice versa. The value of A is determined by the requirement that T c /Td is near 7.8 [15]. This in turn 
determines the value of the constituent mass for all T. 

We expect that the constituent mass m(T = 0) must be substantially larger than the corresponding value for 
fundamental fermions, because a large constituent mass is necessary to delay the onset of chiral symmetry restoration. 
On the other hand, a large constituent mass at the deconfining transitions would be expected to lead to a relatively 
small change in the deconfining temperature. The ratio m(T = 0) /A should be less than one in order for the cutoff 
theory to be meaningful. In the case of fundamental fermions, this ratio is relatively large, on the order of 0.5. We 
have generally found that for adjoint fermions a larger ratio of m(T = 0)/A with T c /Td fixed implies a larger value of 
m(T — 0). We will work with the representative case of mo = and m(T — 0)/A = 0.1. This gives A = 23.22 GeV and 
thus m(T = 0) = 2.322 GeV, with k = 1.2653. For comparison, the critical value of k, k c , below which m(T = 0) = 0, 
IS 7T /8~ 1.234. In figure § we show the constituent mass m and Polyakov expectation value (TtfP) as a function of 
temperature, normalized by dividing by their values at T = and T = oo, respectively. We see that the deconfinement 
temperature Td is very close to its value in the pure gauge theory, due to the large adjoint fermion constituent mass. 
The transition is first order. The constituent mass m has a slow decline to a second-order transition at a substantially 
higher temperature, as indicated by lattice simulations |14l H5] . 



Having established that our PNJL model can provide a useful representation for the finite-temperature behavior of 
SU(3) with both fundamental and adjoint fermions, we turn now to our main interest, adjoint fermions on R 3 x S 1 
with periodic boundary conditions for the fermions in the compact direction. The sole change in the total effective 
potential from the finite-temperature case lies in the contribution of the adjoint fermions, given by Vfl- This single 
change leads to an unusual and unexpected phase structure as L is decreased [2] [9] . 

The effective potential can be written as 



VI. ADJOINT FERMIONS WITH PERIODIC BOUNDARY CONDITIONS 



Veff (m,P) 



V F0 (m) + V FL (m, P) + V g (P) 




(31) 
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Figure 3: The constituent mass m and (TrpP) for two-flavor QCD with adjoint representation fermions with antiperiodic 
boundary conditions as a function of temperature for one choice of k and A. The order parameters are normalized by dividing 
by their values at T = and T = oo, respectively. 



which is general for the case of adjoint fermions with periodic boundary conditions in SU(N). This equation contains 
the basic physics of the phase diagram. The contributions of the first two terms in square brackets are always positive, 
while the final term is always negative. Consider first the n = 1 term 



EH*-** 



{Lrrij) 



M 2 L 2 
2tt 2 



TtaP 
L 4 



(32) 



which is generally the largest term in the series. If either ML is sufficiently large, or some combination of rrijL'a is 
sufficiently small, the overall sign of the expression in square brackets will be positive. The n = 1 term will then favor 
the minimization of TtaP, giving TtaP = — 1 and therefore TtfP = 0. The first case, with L sufficiently large, can 
be identified with the usual low-temperature confined phase of the pure gauge theory. The second case, however, is 
capable of producing a confined phase at small L, provided the constituent mass m is sufficiently small. This will 
certainly occur if chiral symmetry is restored, but it can also occur if L is increased even if m is essentially constant. 
On the other hand, if the sign of the coefficient of TtaP in V e ff is negative, the n = 1 term will favor maximizing 
TtaP, which occurs in the deconfined phase. This simple picture is sufficient for SU(2). The actual phase structure, 
as predicted theoretically and found in lattice simulations, is more complicated when N > 3 [2j|9llIH|- In the case of 
SU(3), there is an additional phase, the skewed phase, in which TrpP is non-zero but negative. As iV increases, the 
number of possible phases increases as well, giving rise to a rich phase structure for small L [9]. 

We consider the behavior of m and TrpP with periodic fermions using the same parameters we used for the 
antiperiodic case. Figure [i] shows the behavior of m and (TrpP) as a function of Lr 1 for the m(L = oo)/A = 0.1 
parameter set, with mo = 0. We see that chiral symmetry breaking persists at L^ 1 = 10 GeV, which is much higher 
than the chiral restoration temperature for antiperiodic fermions. The constituent mass m does fall eventually as L _1 
increases, and chiral symmetry is ultimately restored, but at a temperature on the order of A. In figure[4] TrpP shows 
three distinct phase transitions as a function of L~ l . As L _1 increases, the confined phase gives way to the deconfined 
phase in a first-order phase transition. Because the constituent mass of the fermions is large, the critical value of 
for this transition is approximately equal to T^. As increases, there are two more first-order transitions, from 
the deconfined phase to the skewed phase, and then from the skewed phase to a small-L confined phase we describe 
as reconfined. 

The ordering of the phases seen in the behavior of TrpP for mo = persists as mo is increased. Figure [5] shows 
that with this parameter set the value of L _1 for the confinement-deconfinement transition stays essentially at Tj_ as 
mo increases. The only significant change in the phase diagram is the smooth growth of the extent of the skewed 
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Figure 4: The constituent mass m and (TrpP) for two-flavor QCD with adjoint representation fermions with periodic boundary 
conditions as a function of L _1 for one choice of k and A, with mo = 0. The order parameters are normalized by dividing by 
their values at L = oo and L = 0, respectively. C, D, S and R refer to the confined phase, deconfined phase, skewed phase and 
reconfined phase respectively. 




Figure 5: The phase diagram for two-flavor QCD with adjoint representation fermions with periodic boundary conditions as 
a function of L~ l and Lagrangian mass mo for one choice of k and A. C, D, S and R refer to the confined phase, deconfined 
phase, skewed phase and reconfined phase respectively. The Lagrangian mass is measured in MeV. 



phase as mo increases. Naively, one might expect that the skewed and reconfined phases would disappear as mo — > oo. 
However, at any fixed mo, there will be a sufficiently large value of L _1 such that the fermionic term in the effective 
potential overwhelms the gauge boson contribution, and the skewed and reconfined phases do not disappear, but move 
to very large values of L _1 . On the other hand, taking the limit mo — > oo at fixed L _1 greater than Td always yields 
the deconfined phase. 

PNJL models have a larger set of parameters available, namely {M, mo, k, A}, than in adjoint SU(N) gauge theories, 
which have only mo and A a gj. To the extent that a PNJL model faithfully reproduces the physics of adjoint SU(N), 
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Figure 6: The phase diagram for two-flavor QCD with adjoint representation fermions with periodic boundary conditions in 
the L _1 -k plane for one choice of A. C, D and S refer to the confined phase, deconfined phase, and skewed phase respectively. 

we can think of adjoint SU(N) as giving a two-dimensional surface in the four-dimensional space of PNJL parameters. 
For the case mo = 0, where the Lagrangians are chirally symmetric, we obtain a line through a three-dimensional space. 
This motivates us to consider the phase structure of our PNJL model in the L~ x — n plane with A and M fixed. Note 
however, that the PNJL model will not break chiral symmetry at L^ 1 = unless k > n c — it 2 /8 ~ 1.234. In Figure 
[fj] we show the phase diagram in the L~ x — k plane, with mo = and M and A as before, obtained by numerically 
minimizing V e ff . For most values of k larger than k c , the confined large- L phase and the reconfined phase at small L 
are separated by three phase transitions. There is a transition from the confined phase to the deconfined phase, then 
another transition from the deconfined phase to the skewed phase, followed by a third transition from the skewed phase 
to the reconfined phase. All of these transitions are characterized by abrupt changes in TrpP, while the chiral order 
parameter shows only a slow decrease with increasing temperature. However, there is a narrow range of k between 
approximately 1.250 and k c ~ 1.234 where confinement holds at all temperatures, and chiral symmetry remains 
broken. In this extended phase diagram, the confined and reconfined regions are smoothly connected. Although 
this connection appears only for small range of k values, the corresponding range of constituent mass values is not 
necessarily small. In figure m we show the same phase diagram, now plotted with m(L _1 = 0) replacing k. 

There is a remaining puzzle associated with chiral symmetry restoration. If k < k c , there is no chiral symmetry 
breaking, the adjoint fermions are always light, and TrpP = for all values of The region k < k c , is separated 
from the confined phase with k > k c by a chiral transition in which m > for k > n c . However, for mo light but 
non-zero, there is no true chiral transition across k = k c , only a rapid crossover. The apparent existence within this 
model of a confining phase with unbroken chiral symmetry must give us pause. There is a long-standing argument 
due to Casher that suggests that under very general circumstances confining theories must also break chiral symmetry 
[25] . Confinement without chiral symmetry breaking will also occur in this model if L is sufficiently small, even if 
K > k c - It is not clear if this is due to a defect in Casher's argument, a failure of the PNJL model, or something more 
subtle. 

Figures [6] and [7] show that the small- L and large-L confining phases are connected, at least within the PNJL model. 
A similar smooth connection between the two regions has been observed in lattice simulations using a deformation 
of the pure gauge action in which a TtaP term is added to the action (2]. However, the connection between the 
large-L and small- L regions seen here in the PNJL model may not appear in adjoint QCD without additional terms 
in the action. Our results bear directly on the recent work by Cossu and D'Elia [18 , in which they performed lattice 
simulations of two-flavor SU(3) gauge theory with periodic adjoint fermions. The simulations were carried out on 
16 3 x 4 lattices at various values of the dimensionless lattice parameters, the gauge coupling f3 and the dimensionless 
lattice fermion mass m,;. The size of the compact dimension L in physical units is given by L = 4a, where a is 
determined from (3 and m; via the renormalization group. At fixed mi, larger values of [3 correspond to smaller values 
of L. In these simulations, clear evidence was found for confined, deconfined, skewed, and finally again confined phases 
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Figure 7: The phase diagram for two-flavor QCD with adjoint representation fermions with periodic boundary conditions in 
terms of L _1 and the constituent mass at L" 1 = 0, m(L~ 1 — 0) for one choice of A. C, D and S refer to the confined phase, 
deconfined phase, and skewed phase respectively. The constituent mass is given in MeV. 

as (3 was increased, corresponding to shrinking L. As m/ was decreased, the critical values of /3 at which the three 
transitions occurred moved closer together. However, an extrapolation to m/ = indicates that the two confined 
regions remain separated as m; goes to zero. The deconfined phase clearly persists in this limit, but the skewed phase 
may or may not disappear when mi — 0. All lattice simulations on finite lattices have lower limits on the Lagrangian 
mass mo, so it is a priori possible that the two confining regions seen in the simulations are connected, but for a very 
small range of mo- The phase structure seen in these lattice simulations can also be obtained from our PNJL model 
if the chiral limit, mo = 0, lies near the triple point where the confined, deconfined, and skewed phases meet. For 
example, if we fix m(L _1 = 0) by adjusting k, and then decrease L, we will obtain either one, three, or four phases, 
as may be seen from Fig. [7J The figure shows that the scale over which m(L~ x = 0) must vary in order to obtain the 
different phase structures is quite small, on the order of 150 MeV, compared to the scale of m(L _1 = 0) itself. 

VII. CONCLUSIONS 

We have extended the PNJL treatment of SU(3) gauge theories to the case of adjoint fermions with periodic 
boundary conditions on R 3 x S 1 . This class of models can have Z(N) symmetry for small L provided the constituent 
mass of the fermions is sufficiently light. The constituent mass rather than the Lagrangian mass is relevant because 
the constituent mass determines the fermionic coupling to the Polyakov loop. Previous analyses implicitly assumed 
that the constituent mass could be set arbitrarily, ignoring the effects of chiral symmetry breaking. Because PNJL 
models give a phenomenological description of both deconfinement and chiral symmetry restoration, they are natural 
tools for exploring adjoint fermions with periodic boundary conditions. 

We have shown that a simple model reproduces the known successes of PNJL models for fundamental fermions while 
at the same time reproducing the expected behavior at high temperatures needed with adjoint fermions. The large 
separation between the deconfinement transition and the chiral symmetry restoration transition for adjoint fermion 
theories with antiperiodic boundary conditions requires a PNJL model which reproduces the behavior of the pure 
gauge theory to much smaller values of L than have been considered before. We have extensively studied the case of 
the SU(3) gauge theory with two fermions degenerate in mass. This two-flavor case is particularly simple to work 
with, because the effective potential can be written as a function of a single coupling constant. 

The results for our SU(3) PNJL model with two flavors of periodic adjoint Dirac fermions can be summarized in 
terms of two parameters: k, which represents the dimensionless strength of four-fermion interaction, and mo, the 
Lagrangian mass of the fermions. As with all NJL models, there is a value k c below which chiral symmetry breaking 
does not take place, even as L — > oo. Because chiral symmetry is broken at L _1 = in the corresponding gauge 
theory, k must be taken greater than k c at large L in order to reproduce the behavior of the associated gauge theory. 
For any fixed values of k and mo, there will be a value of L below which the system will be in the confined phase. 
In most of those cases, the system will pass through the confined, deconfined, and skewed phases as L is decreased, 
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before the confined phase is regained. There is a narrow range of parameters where only the deconfined phase lies 
between the two confining regions. If mo is set to zero, there is a small region in the L -1 — k plane, lying above k c , 
that connects the large- L and small- L confined regions. Because the largest contribution to the constituent mass m is 
from chiral symmetry breaking, this behavior will persist for some small range of non-zero mo- Thus there is a single 
confining region, accessible in principle in lattice simulations via a gauge theory to which additional four-fermion 
terms have been added. As we have seen, our results for the PNJL model are completely compatible with the lattice 
simulations of Cossu and D'Elia [18], in which they performed lattice simulations of two-flavor SU(S) gauge theory 
with periodic adjoint fermions. The phase diagram obtained from their simulations can easily be obtained from the 
PNJL model with the appropriate choice of parameters. With more input from lattice simulations of adjoint fermions, 
particularly of basic hadronic properties for large L, we could match the behavior of the PNJL model more closely to 
the gauge theory. 

Adjoint fermions with periodic boundary conditions provide a local field theory that induces confinement in a 
small-L region where semiclassical arguments are valid. It appears likely that the two confinement regions, small-L 
and large-L, can be connected if additional terms, corresponding to irrelevant operators, are added to the action. 
This may have implications for the use of periodic adjoint fermions in achieving volume independence in the large- TV 
limit |26| [27) [28] [29] , which requires unbroken Z( N) symmetry. Perhaps even more important is our newly-gained 
ability to understand confinement in a four-dimensional gauge theory, albeit under somewhat exotic and unexpected 
circumstances. 
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